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Abstract 

In this paper we describe the analytic and numerical evaluation of the 
77 ~^ 1^ process cross section and the Z — )• 777 decay rate within 
the SANC system multi-channel approach at the one-loop accuracy level 
with all masses taking into account. The corresponding package for 
numeric calculations is presented. For checking of the results correct- 
ness we make a comparison with other independent calculations. 
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1 Introduction 



The article describes the implementation into the system SANC [IHS] the scat- 
tering process 

77 7^, (1) 

(see |1H6]) and the decay 

Z 'j'y-f, (2) 

(see [ZllH]) in the Standard Model (SM) at the one- loop level of accuracy 
in i?^-gauge with taking into account of all masses {Z boson and internal 
ones). The processes is interesting from an educational point of view, since 
the calculation of the cross section and the width involves only loops and 
does not contain diagrams of the tree-level and bremsstrahlung, and also 
exploits the multi-channel approach. The work is done in the framework of 
the 4-bosons processes sector extension in the SANC system [9|[T0]. 

In section |5] the multi-channel approach of diagrams calculation is de- 
scribed, when all the particles participating in the process are considered as 
incoming Z777 — )• 0. All the one-loop diagrams, as well as their correspond- 
ing amplitudes in terms of Lorenz expressions (based on the constructed 
basis) and scalar form factors are discussed, as well as the structure of these 
expressions and a proof of zero axial part in the fermionic loop contribu- 
tion to these processes through the application a special sequence of Shouten 
identities. 

In section |3] we discuss the helicity amplitudes, resulting in a chosen 
channel of the process ([1]) or the decay (|2]) as well as formulas for the cross 
section and the decay width, respectively. 

In section H] the package for numerical computations is described. It is 
created on the basis of analytic calculations in the SANC environment. There 
are given the technical instructions and the contents of the package, and 
also the control flags are described. The sanc_4b_vl.00 package contains the 
processes 77 — )■ 77, 77 — )■ 7Z, ([T]) and Z — )■ 777, 02]), see (FigH]), 

In section Owe present the numerical results of the sanc_4b_vl.00 package 
and comparisons with the known in world literature ones. 

In the conclusion we briefly summarize the results of the work. 
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SANC Project Official Home Page Joint InsBtute for Nuclear Research 

SUPPORT OF ANALYTIC AND NUMERIC CALCULATIONS FOR EXPERIMENTS AT COLLIDERS 
SANC Project ~[ New SANC Modules available for download! 

I News I Project Into j Papers I Team Members j Users Support | Download 

SANC Project Download. 

SANC Client - Special interface for networking at SANC Server, 
ii SANC Modules - Stand alone SANC products for precision calculations. 
±] SANC Generators - Stand alone SANC packages for events generation. 

SANC Archives -- archives of previous version products. 

SANC Modules. 
SANC 4b package. 

O 31/07/2012 SANC 4b vl.OO package (105 Kb tgz-file) [last stable version] 

O 19/11/2009 SANC 4A vl.OO package (21 Kb tgz-file) [stable version] 

This package is intended for calculation of the 1-loop full EW correction to four bosons processes like light-by-light scattering, see D. Bardin et al., 
"Light-by-light scattering in SANC", hep-ph/0611ia8. Current processes: 4A scattering, Z3A decay and scattering. 



Copyright © 2002-2008 SANC Project Team. All rights reserved. 
Copyright © Jomt Institute for Nuclear Research (JINR). 

Figure 1: SANC modules download web site [3]. 

2 Precomputation level, channel Z777 — > 

In the SANC system the basic concept of the analytical calculations is precom- 
putation of vacuum building blocks, namely diagrams, in which all external 
particles are considered to be incoming and not lying on the mass shell. 

These are the building blocks that can be used as the elements in the 
calculation of real processes in relevant channels by means of transformation 
of external particles momenta and replacing the squares of momenta by the 
squares of masses. Consider this concept on the example of the process 
Z777 — )■ 0. 

The process at the one-loop level of accuracy is described by two blocks of 
the diagrams with a fermionic and bosonic propagators, respectively. Their 
calculation can be made independently. 

Block of bosonic diagrams consists of three box diagrams, FigJ2] (a), six 
triangular graphs — pinches, Fig|2](b), and three diagrams of the "fish" type 
— self energies, Figj2] (c). 
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Figure 2: Z777 — )• process diagrams. 



Block of fermionic diagrams consists of only three box diagrams, (Figl2] 
(d). Each diagram is characterized by different order of 4-momenta of the 
incoming particle — for Z boson and pi,P2,P3 for photons, respectively. 

In the computation of the diagrams — the application of Feynman rules, 
the Passarino-Veltman (PV) reduction of one-loop integrals [11], scalarizing 
and separation of the poles — the expressions of the amplitude can be rep- 
resented in the form of the sum of products Lorentz structures and relevant 
scalar form factors [12] . 

In terms of Lorentz structures we write the expression for the covariant 
amplitude Z777 — )■ 
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Az,,,^o = 5^ [J-^ (5 , t , n) + J-f (s , t , ] T^^^r (3) 

i=l 

The four rank tensor, with the exclusion by the conservation low of momen- 
tum p4 and imposing the conditions of physical transversality and zero mass 
of the photons {pia = P213 = Psu = and pf = P2 = pI = 0), looks like: 

T"'^'^" = 6af^Pl/3P2u - SapPlfiP2u " Sp^P2aP2u 

+Sl3uPlfiP2a + SfiuP2aP2fi + Sf,uP2aP3l3 

H [5auPl^3Pl^, - 5apPiuP2^, - SafiPl^Plu 

+5auPll3P2fM + 5l3^,PluP2a + S^uPwP3o^ 
+ - {SatiPluPsp - ^auPl^JiP■3P + S[Sf^P2uP3a " ^PuP2tJiP3a - ^/xi/PSaPS/?) 

\ Hi 
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Pl0PlnPluP2a + Pl0PluP2aP2n 

"2 i^a/iPl^Plu + Sa(5PluP2p.) - u5 ^,^PipP2a] - -^Sub afi^ ixv ■, 
t 

PipPllj,P2aP2v - -pipPl^P\vP2a 

+ i^^a^PliiPly - Sal3PlnP2u^ , 
t 

-PlfiPlliPlvP2a + PipP2aP2ij.P2v 

PlfiPluP2aP3fi 

+ {Sa/3PlixP2iy - Saij.Pll3P2iy - (^a//Pl/3P2i/ - bpuP\ixP2a) 
+i ((^a/iPl/3Pli/ - ^a^PlnPlv) " s5a^iPlvPzp\ + 
PluP2aP2fiP3l3 - ^[s {S)3nP2uP3a " 5pvP2^iPZa + SauP2nP3p) 
+tSi3i^PluP2a - U {Sfi^P2aP2v - ^pyP2aP2i,)] + |S^(^ar/(^/3/., 
PllJ.P2aP2uP3l3 " ^ {sSaij.P2uP3P + t5p^Pi^P2a) + ^S^^q^^^^, 
P2aP2nP2uP3P " ^P2^LP2uP3aP30 " (y5l3uP2cxP2n - ^S^^P2f,P3a 
Pll3PlnPluP3a - jP\nPlvP3aP3P " ((^ai^Pl/jPl^i " ^(^at/Pl;iP3^) 
PipPlvP2ij.P3a - ^ {u6auPl^P2n - s5/3/xPli/P3a) + SM^5ay5/3^, 
Pl^PlMP2i/P3a + ^[-S iSauPlixP30 " (^a//Pli/P3/3 " (^/3i.PlMP3a) 
+^ {Sat^PipPlu - SauPl^Pln) - u6a^iPlpP2u] + ^S«^aM^/3i/> 
Pll3P2^iP2uP3a " 2 *<^/3j«^'2i.P3a 
+ ((^a/3Pli/P2/^ - SauPll3P2n - ^PnP\uP2a) 
+ {Sf3,j.P2aP2u - Sa(3P2,,P2u) + ^StSauS/S,,, 



Tii^" = Pl^,P2uP3aP3l3 " ^ {t^ l3uPltiP3a + u5a^,P2uP3p) + ^tU^af^S/s^. 

The precomputation files AAAZ Box, AAAZ pinch, AAAZ fish (see the 
SANC process tree in Figj3] from SANC chent [3]) contains the sequence of 
procedures for calculation of the covariant amplitude. 

Form factors J-^ are the scalar coefficients in front of basis structures of 
the covariant amplitude. They are presented as some combinations of scalar 
PV functions Aq, Bq, Cq, Dq yjj, and depend on invariants s ,t ,u , and also 
on fermion and boson masses. They do not contain ultraviolet poles. The 
derived one-loop scalar form factors can be used for any cross channel after 
an appropriate permutation of their arguments s, t, u. 

Explicit expressions for the boson and fermion parts of the forms factors 
are not shown in this article because they are very cumbersome. A complete 
answer for J-'j one can be found in the package which is downloadable from 
the homepages of the computer system SANC. Note that the expression for 
the amplitude of boson diagrams are similar to fermion one except for the 
explicit representation of form factors. 

To the Lorentz structure of the expression, the axial interaction of Z- 
boson with fermions g^f^ai3ufi does not give the contribution due to charge 
symmetry. One can be shown that it cancels in a full set of diagrams. 

For the analytical proof of this fact a special sequence of the Shouten 
identities was used: 

where /Xj? denotes any index. 

Namely, we applied, step-by-step, the basic identity |5] contracted with 
certain number of 4-momenta. It was found only 5 families of the Ihs of the 
Shouten identities for the substitutions. 

• Family epj^p^p^aiSfi^ifia? has 4 members: 

Cpip2P3"'^M5?,M6?' ^PlP2P3l3^afJ.5?-i ^piP2P3/3'^'^M6?' ^p\p2P3il^U^l(,l ■ 

• Family ep^7^p^?^al,j3l5^J.'!,v'! has 3 members: 



6 



• Family e^. a?,i/?Pfe/i? has 11 members with 3 sets each 

(first set for i,j — 1,2, second set for i,j — 1,3 and the last one for 

i,i = 3,2): 

^PiPjoivPZfii ^p^pja'!fiP3l3i ^piPj/3?uP3ai 

^PiPjauPltj,P3l3, ^piPjafiPluP3l3, ^piPjPnP2vPZa-, ^piPjl3uP2fiP3a, 
^PiPjl3i>7P2aj ^piPjO.v'fPllS-i ^piPjixvPlfSi ^piPjl3vlP2a-i 

• Family name ep-??7Pk^? has 3 members with 3 sets (for z = 1, 2, 3): 

^Pi^HvPSat ^piP^,uP2on ^pianvPl/S ■ 

The procedure Pro ject ions () was call in between. 

#procedure Projections () 
id pl(al)=0; 
id p2(be)=0; 
id p3(nu)=0; 

id p3 (mu) =-pl (mu) -p2 (mu) ; 

repeat id pl.pl=0; 
repeat id p2.p2=0; 
repeat id p3.p3=0; 
repeat id pl.p2=-s/2; 
repeat id pl.p3=-u/2; 
repeat id p2.p3=-t/2; 
#endprocedure 

3 Processes level, helicity amplitude. 

When we implement any processes, we create the building block for annihila- 
tion to the vacuum and then use these building blocks several times replacing 
incoming momenta p's by corresponding kinematic momenta with the right 
signs. 

The covariant amplitude for the channel 77 — )■ 7Z can be obtained from 
annihilation to the vacuum with the following permutation of the 4-momenta: 

Pi Pi, (5) 

P2 P2, 

P3 -P3, 
Pi -P4 



7 



and for the decay Z — )■ 777 are: 



P4 P4- 



(6) 



Furthermore we calculate J-^ by module AA — i- AZ (FF) , Z — t- AAA (FF) , 
then helicity amplitudes by the module AA — )■ AZ (HA) , Z — )■ AAA (HA) and 
finally — the analytic expression for differential and total process cross sec- 
tion in sanc_4b_vl.00 package. 

^ Root 
¥ SANC 

0- □ OED 

■a EW 

'\ 'C3 Precomputation 
¥ Q Self 

¥ Q Boson 

9- □ Ren Blocks 
^ Boson Self 
Tadpole 
- AAAA fisfi 
AAAZ fisfi 
□ Fermion 
3 Vertex 
o- □ bbb 
o- □ bff 
9 Q bbbb 

|- AAAA pincfi 
\w AAAZ pinch 
t "3 Box 
0- □ ffff 
0- □ ffbb 
9 Q bbbb 

AAAA Box 
'-fe AAAZ Box 
□ Processes 
Ifi QCD If □ QCD 



? |(t' SANC 
5>- □ QED 
t Q EW 

f □ Precomputation 
f 3 Processes 
«- □ 3 legs 
t 3 4 legs 
«- □ 4f 
o- □ 2f2b 
t Q 4b 

I) Q AA -> AA 



BlAA-> AA(FF) 
ItoAA-> AA(HA) 
a AA -> AZ 
I fTOAA-> AZ (FF) 
L ^ AA -> AZ (HA) 
Z -> AAA 
fwZ-> AAA (FF) 
BlZ -> AAA (HA) 



Figure 3: Z777 — > SANC processes tree. 



For more compact presentation of the results and more effective numeri- 
cal implementation in the system SANC it was applied the method of helicity 



8 



amplitudes. Lorentz structure of the classical expression is contracted with 
polarization vectors, and one gets the orthogonal set of scalar variables, ex- 
pressed in terms of forms factors — helicity amplitudes [T3] . 

The expressions for helicity amplitude with substitution of boson and 
fermion forms factors for both channels were calculated. It was received the 
full agreement as compared with the calculated in the literature (though 
some trivial misprints had to be corrected). 

The cross section of the reaction 77 — )■ 7Z is calculated by the formulae: 

da^-y^-yZ ~ , = |^77^>7Z| d^^ \ 

where A^^^^z is the covariant amplitude of the process and d^^"^^ is the 
two-body phase space: 

' = {2tt) 6 [pi +P2-P3- Pa) 



(27r)3 (27r)3 ' 
For the differential cross section one gets: 

da^^^^z = — — ( 1 ) \A^y^^zf d cos 6, 

where 6 is the scattering angle of the Z boson in the center of mass system. 
The decay width of the Z boson is calculated by the formula: 



3!384^^3M| J \^z^^^^\^dsdtdux6{Ml-s-t 



u 



4 Technical description of the package 

The process ([1]) and the decay ([2]) and also previously implemented process 
77 — )■ 77 are included into the package sanc_4b_vl.00 (Fig|3]), which can be 
downloaded from the project site [3]. 

Our main purpose is to introduce the package sanc_4b_vl.00 into generic 
integrator SANG, based on an algorithm VEGAS [TJJ. 

Here we present the technical description of this package — main flags 
and the options: 

1. bbbb_main.F — the main file. 
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2. bbbb_ha_ll_ll.F — the HA of 77 — > 77 process from SANC system, 

3. bbbb_ha_ll_12.F — the HA of 77 7Z process from SANC system, 

4. bbbb_ha_2_lll.F — the HA of Z — >■ 777 process from SANC system, 

5. *.f — the hbrary of special functions and algorithms, 

6. *_input.h — the set of various setups of input parameters, 

7. README, INSTALL and other instructions files. 

In README and INSTALL files one can find instructions how to use the 
package. The main options one can change in bbbb_main.F. 

pid(I) — choice of the process: 

• I = AA2AA, 77 — )■ 77 process 

• I = AA2AZ, 77 7Z process 

• I = Z2AAA, Z 777 decay 

ipm(I) — choice of incoming photons helicities sum in cross section: 

• I = SS, total helicities sum 

• I = ++, " + +" helicities sum 

• I = H — , " H — " hehcities sum 

itl(I) — choice of helicities for Z-boson: 

• I = T, Z-transverse helicities 

• I = L, Z-longitudinal helicities 

• I = S, Z-helicities sum 

iqed(I) — choice of calculations for QED (fermionic) corrections: 

• 1 = 0, without QED corrections 

• 1=1, with QED corrections 
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iew(I) — choice of calculations for EW (bosonic) corrections: 

• 1 = 0, without EW corrections 

• 1=1, with EW corrections 

To get full EW answer with the interference one should set iqed = 1 , iew 
= 1. 

gfscheme(I) — choice of the EW scheme: 

• 1 = 0, ao calculation scheme 

• 1=1, Gp scheme 

• 1 = 2, G'f scheme, 

when ao is replaced by aoj, = V^GpM^ (1 - M^/M|) /tt. 

isetup(I) — choice of the setup: 

• 1 = 0, Standard SANC input [PDG 2006] 

• 1=1, Les Houches Workshop (2005) 

• 1 = 2, Tevatron-for-LHC Workshop (2006) 

• 1 = 3, Custom setup 

5 Numerical results 

The calculations were made for the cross-section of the reaction 77 — ?> 7Z 
and decay width of the Z 777 with the following values of parameters: 

q; = 1/128; 
TT/6<e < 57r/6; 
Mw = 80.22 GeV; 
Mz = 91.173 GeV; 

me = 0.1 GeV, = 0.1 GeV, = 0.1 GeV; 

m„ = 0.1GcV, mc = 0.1GeV, = 120.0 GcV; 

ma = 0.1 GeV, = 0.1 GeV, rub = 5.0 GeV. (7) 
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The numerical results for the reaction 77 — )■ 7Z were compared to [IHH] 
componentwise for boson, fermion contributions and their interference taking 
into account the helicities of Z boson, and initial photons (for ++ see FigJU 
for H — see FigjS]) in the energy range from 100 GeV up to 2 TeV. 



++ incoming photons helicities 




200 400 600 800 1000 1200 1400 1600 1800 2000 
Vs, GeV 



Figure 4: 77 — )■ 7Z SANC cross section (" + +") incoming photons helicities. 



+- incoming photons helicities 



10"^ 









1 fermions loop Z longitude 




bosons loop Z longitude 




p--.,..,^^ interference sum Z longitude 




f ""---^^^^Trt-^^.^ fermions loop Z transverse _ 




""■^^^^s;,..,.^^ bosons loop Z transverse 




^=^^>s,,^^ interference sum Z transverse 







200 400 600 800 1000 1200 1400 1600 1800 2000 
Vs, GeV 



Figure 5: 77 — )■ 7Z SANC cross section (" H — ") incoming photons helicities. 

For all of the contributions it was obtained a good agreement with the 
results, given in the literature. 
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When calculating the decay width Tz, the variation of the cut parameters 
(the angle 9 cut and the energy of photons A/i^nt) performed in a large 
range of values (Fig|6]). It was found a wide plateau of stability, giving a 
result consistent with those given in the literature 0. 




Figure 6: Z — )■ 777 SANC decay width plateau of stability. 



6 Conclusions 

The precomputation strategy of the SANC system ^ and the place of the 
process ([1]) and ([2]) on the SANC process tree were presented. 

The implementation of analytical results and the concept of modules were 
described. 

Its numerical results of the sanc_4b_vl.00 package were compared with 
those existing in the literature. The package is available for download at web 
page [Sj. 

The authors are grateful to S. Bondarenko, G. Nanava, V. Kolesnikov, 
and A. Sapronov for useful discussions of numerical and analytical calcula- 
tions. 
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